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Abstract. We investigate a circuit cavity electromechanical system, where two 
superconducting coplanar waveguide cavities (CPWC's) are capacitively coupled by 
a nanoscale mechanical resonator (MR). We show that, with this electromechanical 
system, two-mode continuous-variable entangled states of microwave electromagnetic 
fields can be synthesized deterministically either via coherent control on the dynamics 
of the system, or through a dissipative quantum dynamical process. The first scheme 
operates in the strong coupling regime, and is insensitive to the initial thermal 
noise, while the second one exploits the mechanical dissipation as a useful resource 
and only needs high frequency low-Q mechanical oscillators. These protocols may 
open up promising perspectives for quantum information processing with cavity 
electromechanics. 
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1. Introduction 

Circuit cavity electromechanics [H El El H], the counterpart of cavity optomechanics 
[51 [6] in the form of electrical circuits, describes the parametrical coupling between the 
motion of a micro or nanoscale MR and an electrical circuit. The underlying physics 
of cavity electromechanics is that the motion of the mechanical oscillator modulates 
the capacitance of the electrical circuit, thus creating parametrical coupling between 
these two systems. Compared to their optomechanical analogues, electromechanical 
systems have the advantages that these low-loss superconducting circuits are easily 
cooled to ultralow temperatures [7], and can be fabricated on a single chip using 
the standard optical lithographic techniques. Recent experimental and theoretical 
progress has shown that cavity optomechanics and electromechanics are pretty useful for 
macroscopic tests of the fundamental laws of quantum mechanics, or for other practical 
applications relevant with quantum phenomena [H [91 [101 ECU El2l 1131 [HI HSl [161 CEO 
UHl H9l ETil EH [221 [231 EH [251 EH [281 EH1 EOl EI]. Of particular interest is the 
generation of non-classical motional, photonic and hybrid quantum states for basic 
tests of quantum theory, as well as applications in quantum information processing 

[321 E31 EH EHl [361 EZ1 EEl EHl SQl SH S21 S31 SHI ffl ST]. 

In this work, we study a circuit cavity electromechanical system consisting of 
two superconducting CPWC's and a nanoscale or micro MR. In this proposal, the 
superconducting cavities are capacitively coupled by a capacitor that incorporates the 
nanoscale or micro MR into its electrode plates, and is biased by a driving voltage. 
We show that, through suitably choosing the driving frequencies of the voltages, we 
can generate various linear operations between cavity photons and mechanical phonons 
on demand. In particular, with this circuit electromechanical system we present two 
different protocols for synthesizing continuous-variable entangled states of microwave 
electromagnetic fields based on currently available experimental setups. The first 
protocol, operating in the strong coupling regime, relies on coherent control over the 
dynamics of the system and is insensitive to the initial thermal noise. The second one 
is based on a dissipative quantum dynamical process, which exploits the mechanical 
dissipation as a resource and only needs high frequency low-Q mechanical oscillators. 
These protocols may have promising applications for continuous-variable quantum 
information processing with cavity electromechanics. 

2. Coupling two superconducting CPWC's via a MR 

As shown in Fig.l, we consider the electromechanical system where two superconducting 
CPWC's are capacitively coupled by a nanoscale or micro MR, which is driven by a gate 
voltage V x . To implement this scheme, the superconducting CPWC's are electrically 
connected with a capacitor biased by a driving voltage. The capacitor is formed by 
two parallel metal plates, one of which is replaced by a metallic membrane realizing 
the MR (drum resonator) [2]. Alternatively, one can choose to use a micromechanical 
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Figure 1. (Color online) The schematic of two superconducting CPWC's capacitively 
coupled by a nanoscale MR, which is driven by a gate voltage V x . 

bulk dilatational resonator, as recently used in the experiment to couple with a phase 
qubit [7]. The working frequency of these MR's is in the range of GHz, and they can 
couple to the CPWC's through interdigitated capacitors. In both cases, the capacitance 
C of the capacitor is dependent on the MR displacement X = \JWj (2muj m )(b + tf), 
where m is the mass of the MR, uj m the mechanical vibration frequency, and b the 
annihilation operator for the MR. If we assume that the displacement is much smaller 
than the equilibrium distance d between the metallic membrane and the metallic base 
electrode, then the capacitance approximately becomes C = Co(l + X/d), where Co is 
the capacitance for the MR in equilibrium. 

For a superconducting CPWC [IE], the voltage at position x is 

Vj(x) = J^(a] + d 3 ) cos(27rx/L,), (j = 1,2), (1) 

where Uj is the resonant frequency, Cj the total capacitance, dj the annihilation operator, 
and Lj the length for the jth CPWC respectively. With a coupling capacitance C 
between these CPWC's, the coupled interaction can be derived as 

Ht = ^C (l + X/d){VM + V 2 (0) - V x f. (2) 

We subsequently perform a rotating-wave approximation to simplify the coupled 
interaction. After neglecting rapidly oscillating and other higher order terms, the 
Hamiltonian describing the coupled system can be derived as 

M 3 = hLOid\di + huj2a\a2 + hio m b^b 

- h gi (t)(b + 6 t )(ai + a\) - hg 2 (t)(b + tf)(a 2 + a J) (3) 

where gj(t) = ^ \j 2~m^f~tT^ x (*) ' anc ^ we have included the free Hamiltonian of the two 
CPWC's and MR in the first three terms. The last two terms in Eq. ([3]) are the 
linear interaction between the CPWC's and the MR. Up to now the result is valid 
for arbitrary driving voltage signals. Through adjusting the driving frequency of the 
voltage, we can generate various linear operations. For instance, the beam-splitter 
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interaction between two cavity modes and the mechanical mode, as discussed in Ref. 
[49. 50 i to realize intracavity state transfer, can be recovered from Hamiltonian (3) 
by choosing V x (t) = V^coscOdt and setting u>d = u\ — u m = w 2 — u m . In this way, 
the modulation of the coupling strength provides an effective tool for controlling the 
interaction between the CPWC's and the MR. In the following section, we will show 
how to prepare two-mode squeezed states of the CPWC's through engineering the desired 
interaction between photons and phonons. 



3. Generating continuous- variable entanglement between the CPWC's 

We now consider the case where the MR is driven by a gate voltage of the form 
V x {t) = V x l cos(a^t) + V^ 2 cos(cu^ + (ft), where is a fixed phase difference between the 
voltage components. If we choose the driving frequencies as uj\ = ui+u m , uj\ = \co m —uj2\, 
corresponding to the blue sideband and red sideband driving for the MR, then under the 
rotating-wave approximation we can obtain the Hamiltonian in the interaction picture 

= -ne^alP + oiS) - ne 2 (alb + d 2 P), (4) 

where 

The Hamiltonian (jlj) describes a system of three coupled harmonic oscillators with 
controllable coefficients [511 [52j [53]. The first term describes simultaneous creation or 
annihilation of a photon in CPWC1 and a phonon and is responsible for entangling the 
CPWC1 and the MR, while the second term describes the exchange of excitation quanta 
between the CPWC2 and the motion. These terms together will lead the CPWC's to 
be entangled with each other. In what follows, we will discuss two different schemes 
to realize this goal, one of which is based on coherent control on the dynamics of the 
system, while the other is through a dissipative quantum dynamical process. 



3.1. Dynamical generation of two -mode squeezed states 

We first focus on the regime where the dissipative effects on the coherent dynamics can 
be neglected, i.e., the strong coupling regime, {9i, 9 2 } 3> {k>i, K2i n thlm-,lm\i where Kj 
is the jth CPWC field decay rate, nth is the thermal equilibrium occupation number 
for the mechanical mode at temperature T, and 7 m is mechanical dissipation rate. 
This regime can be easily realized, since the coupling strength Qj can be tuned by the 
classical driving amplitude Vj , and high-Q superconducting CPWC's and MR's can be 
conveniently fabricated in the laboratory. In this limit the coherent dynamics of the 
coupled system can be easily solved in the Heisenberg representation. 
The Heisenberg equations of motion read 

ai = zeifr, (5) 

a 2 = iQ 2 b, (6) 
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b = i®ia\ + i6 2 a 2 , (7) 

which would generate periodic dynamics provided that | O2 1 > After some 

straightforward derivations, we can obtain the time evolution of the operators as 

Ht) = ^St(O) sinBt + ^[1 - cos6t]4(0) + ^[|6 2 | 2 

- |6i| 2 coset]ai(0), (8) 

a *(*) = sin0t - ^rl 1 - coset]oj(0) - ^[IGil 2 

- |6 2 | 2 coset]a 2 (0), (9) 
b(t) = 6(0)cos6t+ ^[i6 2 a 2 (0) +i6ia t 1 (0)]sinet, (10) 



with 6 = y |@ 2 | 2 — |@i| 2 - In general these solutions describe tripartite entanglement 
among cavity modes and the mechanical oscillator. However, we find at the instant 
T n = tt/Q Eqs. (|8])-(II0|) become 



|0il 2 + l© 2 | 2 . , n . , 26x62 

-Q2 a l(°) + -@2 



ai(^r) = 7^ + -7^ a 2(0)> (11) 



a 2 (T„) = - I9l|2 + I92|2 a 2 (0) - ^at(O), (12) 
b(T n ) = -6(0). (13) 

Therefore, at this instant the motion is decoupled from the cavity modes and returns 
to its initial state. Moveover, at the time T n the two cavity modes are entangled with 
each other. 

To be more specific, we introduce the unitary operator S(() = e C i a 2-Caia 2) where 
£ = tanh~ [2r/(l + r 2 )],r = |6 2 /6i|. Then we find Eqs. (1111) and (jl~2l) can be rewritten 

as 

ai(T w ) = coshCai(O) + sinhC4(0) = Sai(0)&, (14) 
a 2 {T n ) = - (coshCa 2 (0) + sinhC4( )) = e in Sa 2 (0)Sl (15) 

These results imply that at the instant T n the time evolution operator of the whole 
system corresponds to U(T n ) = S(CV ® I m , where I m is the identity operator for 
the mechanical mode. Thus, if initially the MR density matrix is a thermal state at 
temperature T given by 

g m (0) = (1 - e -^/k B Ty-H m /k B T ( 16 ) 

where ks is the Boltzmann constant, and H m = fkJmfob + 1/2), then Q m iT % ) = g m (Q) 
in the Schrodinger representation. Moveover, at the time T n , the two cavity modes will 
be prepared in a two-mode squeezed state if the initial state is the vacuum state for 
both cavity modes, |00) c . In particular, using the factored form of the two- mode squeeze 
operator [51] 

S(C) = ( cos } 1 ^ _1 e _a i a 2 tanh f e ~( a i a l+4 a 2) ln ( cosh f) e a i a 2tanhC fyj\ 
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at the time T w the state of the two cavity modes is 

^ n=0 

= (tt^)E(tt^J i"' B >- (18) 

x ' n=0 

The state (fl8|) is a two-mode squeezed state of the photon fields in the two cavities, which 
exhibits Einstein-Podolsky- Rosen (EPR) entanglement. The mechanical motion of the 
MR plays a fundamental role in establishing the entanglement, nevertheless the initial 
motional state does not affect the efficiency of the scheme. The degree of squeezing 
(squeeze parameter () is determined by the ratio of G2 to Gi, which can be controlled 
on demand through tuning the driving signals. When the squeezed state is generated 
at the time of T n , we switch off the couplings between the CPWC's and the MR. Then 
the squeezed state can be preserved until the cavity fields are coupled out. 




Figure 2. (Color online) The total variance V versus 8t for different values of the 
parameters n th and r: (a) n th = 0, 10, 10 2 , r = 1 + 2.5~ 5 ; (b) r = 1 + 2.5~ 2 , 1 + 
2.5~ 3 ,l + 2.5- 5 , n t h=10. 



In order to check the above analysis, we exploit the total variance V = ((Am) 2 + 
(Av) 2 ) of a pair of EPR-like operators u = X\ — X 2 , and v = J\ + Pi [55], with 
Xj = {cij + at) /a/2, and Pj = —i(a,j — a])/\/2,j = 1,2. According to Ref. [55], a 
two-mode Gaussian state is entangled if and only if V < 2. For the two-mode squeezed 
vacuum state S' t (C)|00) c , the total variance V = ((Am) 2 + (A-0) 2 ) = 2e 2( % implying 
this state exhibits EPR entanglement. In Fig. 2, the quantity V is plotted versus the 
scaled time 0t, for different values of the initial thermal phonon number n t h and the 
parameter r. Fig. 2(a) shows the total variance V as a function of the scaled time Qt 
under different values of the thermal phonon number n t h for a fixed parameter r. From 
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Fig. 2(a) it can be found that the total variance V is unaffected by thermal noise at the 
instant T n , i.e., at half period it is independent on the initial thermal phonon number 
nth- This result can be understood since our discussions hold under the conditions that 
the couplings of the system to the thermal reservoirs can be neglected. In such a case 
the coherent dynamics is purely governed by the Hamiltonian flU) and thermal effects 
enter only as an initial condition for the mechanical mode. Fig. 2(b) displays V versus 
Qt under different values of the ratio r with a fixed mean phonon number nth- We find 
that at half period the two CPWC's are steered into a two-mode squeezed state starting 
from a thermal state for the MR. The degree of squeezing and the total variance for the 
two-mode squeezed state is determined by the parameter r. These numerical results are 
in accordance with the analytical conclusions above. 

3.2. Two-mode squeezed states at steady state through a dissipative dynamical process 

In the previous section, we have discussed how to prepare a two-mode squeezed 
state for the CPWC's via coherent control on the evolution of the coupled system. 
Though the protocol seems promising, the experimental implementation of this scheme 
requires stringent conditions, i.e., the coupling to the environment reservoirs should be 
neglected, thus requiring very high quality factors for both the CPWC's and MR's. For 
superconducting CPWC's, this condition can be fulfilled since high-Q superconducting 
stripline cavities are easy to be fabricated in the laboratory. However, as for the 
mechanical oscillator, in particular which has to be incorporated into a capacitor in 
order to couple to the CPWC's, this requirement is too demanding. In addition, high 
frequency MR's are required when the quantum regime is entered. However at present 
these GHz mechanical oscillators are plagued by very low quality factors. It is known 
that MR performance degrades considerably as the oscillating frequency increases. 
In this section, we will present an alternative scheme which exploits the mechanical 
dissipation as a useful resource and only needs high frequency low-Q mechanical 
oscillators. 

In what follows, the system-environment interaction is assumed Markovian, and 
then is described by a master equation in Lindblad form. We assume that the CPWC's 
couple with the vacuum bath, but the MR couples with a thermal bath. Then the time 
evolution of the density operator p for the whole system is described by the master 
equation 



dp 
~dt 



[je, p\ + jSf Cl/ 5 + J*f C2 P + Jf m p, 



(19) 



where 



(2djpOj — djdjp — pdjdj), 




+ ^n th (2ppb - bPp - pbb ] ). 



(20) 
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In the following we focus on the regime where 7 m 3> {0i, 62} 3> {acx, k 2 , nthjm}- The 
condition 7 m ^> {©1,02} corresponds to strong mechanical damping for the MR, i.e., 
very low quality factors, while ntnlm 7m implies near zero temperature for the 
mechanical modes, thus requiring ground state cooling of the MR. In effect, in the 
regime of large mechanical frequency (in the GHz range) and at cryogenic temperature, 
the thermal phonon number is nearly zero, i.e., n t h = ( e riw rn/k B T _ -q-i ^ w j 1 j c j 1 
corresponds to coupling with the vacuum bath for the MR. Under this regime the 
master equation ( |T9i) then can be approximated as 

f t = - l[je,p) + ^{2b P P - Pb P - P b% (21) 

We now introduce the phonon number representation for the density operator p, 
i.e., p = ^ n=0 p mn |m)(n|, where p mn are the density-matrix elements in the basis of 
the phonon number states {\n),n = 0, 1,2, ...}, and are still operators with respect to 
the cavity fields. Under the condition of strong mechanical damping, the populations 
of the highly excited motional states can be neglected. Therefore, we consider only the 
matrix elements p mn inside the subspace {|0), |1)} of the phonon numbers. In this case, 
the master equation (I2TT) leads to the following set of coupled equations of motion for 
the density-matrix elements 

Poo = - 2poi(@iaI + ©2^2) + i(Qiai + 2 4)Pio + TmPu, (22) 

P11 = - «Pio(©iOi + + i(Qia{ + 6 2 a 2 )poi - JmPn, (23) 

P01 = - «Poo(@iai + ©24) + i(Biai + @2a 2 )Pn - ~yA)i- (24) 

In the regime of strong damping rate j m , the elements poi an d P11 can be adiabatically 
eliminated from the above equations [56], leading to 

p aL =—(&Pii-p 00 &) (25) 

7m 

with 

9 = ^ + ( 26 ) 
The reduced density operator for the CPWC's can be approximated as g c = Tr m (p) ~ 
Poo + Pn- Replacing ( 1251) into ( 1221) and ( 123]) . and adding up them, after neglecting 
higher-order terms we obtain the evolution of the cavity modes with an effective master 
equation 

^ = ^ (2&Q C & - &®Q C - Q C &@) (27) 

with T c = 40 2 /7 m . This master equation has the form of the standard engineering 
reservoir scheme. The only pure steady state of the system is the eigenstate \ip) of the 
operator @ with zero eigenvalue, ensuring that there is no further eigenstate \<f>) of @ 
such that [@,<&fi]\<f>) = 0. For the operator @ = ^a 2 + %-a{, we will find that this 
condition cannot be satisfied. From the eigenvalue equation 3>\ip) = 0, and the relation 
Q = Sd 2 S^, with S(s) = e 5 * 1 * 3- ^, ? = tanh _1 [ei/e 2 ], one can readily find that 

\^) e = S\fi,0) c (28) 



Synthesizing continuous-variable entangled states of microwave fields with cavity electromechanics9 



is a steady state of the master equation (1271) . but not the only one [57J EB]- Here \p, 0) c 
denotes an arbitrary state for the first CPWC mode and the vacuum state for the second 
one. 

In order to steer the system into the two-mode squeezed state <S|0,0) C , we need 
another dissipative process together with the described one, leading to the effective 
master equation 



with the operator 3 = Sa,iS\ This master equation describes simultaneous ground 
state cooling of the system in the transformed picture with the basis 3, 3 [59, 60J. In 
fact, one can find that #5|0,0) c = ^<S|0,0) C = 0, and [3, 3} = 0. Therefore, the 
unique steady state of master equation ( 129]) is just <S|0, 0) c . The Lindblad term with 
respect to the 3 operator can be engineered from the Hamiltonian 



following the same reasoning as that for the 3 operator, with the driving frequencies 
chosen as uj\ = \u m — coi\, and cu% = co 2 + oj m . However, for the case of just one MR, 
we cannot get the master equation ( 129]) to realize simultaneous cooling of both modes 
3 and 3, since both and M" cannot be possessed simultaneously only through 
adjusting the driving frequencies for one MR. In order to have both cooling processes, 
in 3 and 3, we can employ a stroboscopic cooling scheme. In this approach, the system 
evolves during a time t in N cycles of duration St = t/N, while the driving parameters 
alternate between the ones with respect to 3 and those of 3. The stroboscopic limit 
is valid provided that the time interval 5t is much smaller than 1/r, in which case the 
effective dynamics of the system is just as that described by the master equation ( 129]) . 
Alternatively, one can couple the CPWC's with two MR's, each of which is driven by 
a bichromatic microwave signal to induce sidebands in the CPWC-MR coupling. In 
this case, one can realize the Hamiltonian M' for one MR, and simultaneously have 
the Hamiltonian M" for the other. In the regime of strong mechanical damping for 
both MR's, one can exploit the engineering reservoir scheme to get the effective master 
equation ( 129]) . 

It is necessary to verify the model through numerical simulations. To provide an 
example, here we consider the two-MR case, where the dynamics of the system can be 
simulated by the following master equation 



dg c T, 
~dl ~ ~2 




(29) 



(30) 



dp 



\3tf, p] + £> Cl p + Jz? C2 p + ^ mi P + 3?mJ, 



(31) 



with 



— hQ 1 (d\b\ + ai&i) — hQ 2 (°V ) i + ^2^1) 

- frOi(a|6 2 + olibl) - hB 2 (d 2 bl + d 2 b 2 ) 



(32) 



^(n th +l)(2b j p^-^b j p-p^b j ) 
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+ J^nth&lpbj - b 3 b]p - pb 3 b]). (33) 




1 1 1 1 1 1 ' 

12 3 4 

time (1/r c ) 

Figure 3. (Color online) Plot of time evolution of the total variance V through 
numerically solving the master equation pip , together with the result for an ideal 
two-mode squeezed vacuum state. The relevant parameters are chosen as 02 = 2©i, 
7mi - 7m 2 — 7m = 150i, n th = 0.01. 

In Fig. 3 we illustrate the time evolution of the total variance V together with the 
result for an ideal two-mode squeezed vacuum state, under different values for the decay 
rate k, where K\ ~ «2 — « is assumed. The initial state of the system is chosen as the 
ground states for the cavity modes and the mechanical modes. The related parameters 
are chosen in such a way that they are within the parameter range for which this scheme 
is valid and are accessible with current experimental setups. From this figure we find 
that, when the effective decay rate for the engineered reservoir is much larger than 
the decay rate of the cavity photons, i.e., T c 3> k, at steady state nearly ideal EPR 
entanglement (O2/Q1 — 2,V — 0.667) between the photon fields in the two CPWC's 
can be established with a fidelity of 99.98%. The time for reaching the stationary state 
is about T ~ 4/T c . However, when k is of the same order of magnitude of or larger than 
r c , the ideal EPR entanglement is severely spoiled, since in this case the influence of 
the natural reservoir on the generated two-mode squeezed state is almost the same as 
that of the engineered reservoir. Thus, to reduce the influence of the natural reservoir 
requires T c ^> k, or 40 2 ^> 7 m re. To make sure the protocol is valid, we also require 
that 7 m 3> 0. These conditions together imply that 7 m >9> -y/7m K /2. 

4. Implementation 

Regarding the experimental feasibility of the proposals, currently available experimental 
setups of cavity electromechanics [2j [8] are promising platforms for realizing the schemes. 
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We consider superconducting CPWC's with the fundamental frequency of 2n x 10 GHz, 
whose damping rate can be as low as k/27t ~ 10 kHz given a quality factor Q = 10 6 from 
recent circuit QED experiments. The vacuum-fluctuations-induced voltage between the 
central conductor and the ground plane of the CPWC's is typically of order of /iV. As for 
the nanoscale MR's, we can choose to utilize an aluminium membrane integrated into a 
capacitor [2] for the first scheme, or a piezoelectric dilatation resonator [7J for the second 
one, which comprises a piezoelectric thin film of aluminium nitride, sandwiched between 
two aluminium metal electrodes. The superconducting CPWC's and nano MR's can be 
fabricated on a single chip with wafer-scale optical lithographic techniques. For a nearly 
circular membrane with a diameter of 15 /im and a thickness of 100 nm [2], drum-like 
modes are allowed to resonate freely. The fundamental mode is u m /27T = 10.69 MHz, 
giving a zero-point motion of 4.1 fm and a damping rate , y rn /2iT = 30 Hz pj. For the 
first scheme, with the chosen parameters V* = 10 V, V% = l.OlV^, d = 50 nm, Cq = 40 
fF, we get ~ 3 MHz, and the operation time for generating the target state is about 
T n ~ 1 fis. This time is much shorter than the photon life time, and the decoherence 
time for the mechanical mode with about 10 3 phonons. With regard to the second 
scheme, the piezoelectric dilatation resonator is particularly suitable [TJ, which has 6 
GHz frequency and very strong damping rate 'y m /2ir ~ 23 MHz (Q ~ 260). At the 
temperature T ~ 25 mK, the number of thermal phonons in the mechanical mode is 
less than 0.07. If we assume that V* = 1 V, V,? = 2V*, d = 50 nm, C = 25 fF, then we 
have O ~ 17 MHz. The time for reaching the stationary state is about 4/T c ~ 0.5 fis. 

5. Conclusions 

To conclude, we have investigated a circuit cavity electromechanical system, in which 
two CPWC's are capacitively coupled by a MR. With this cavity electromechanical 
system, we have presented two different schemes to generate two-mode continuous- 
variable entangled states of microwave electromagnetic fields confined in the cavities. 
The first scheme, operating in the strong coupling regime, is based on coherent control 
over the dynamics of the system, and is insensitive to the initial thermal noise. The 
second one is based on a dissipative quantum dynamical process, which exploits the 
mechanical dissipation as a useful resource and only needs high frequency low-Q 
mechanical oscillators. These protocols may have interesting applications in quantum 
information processing with electromechanical systems. 
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